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REAL-TIME DIGITAL-COMPUTER-HARDWARE 

SIMULATION OF A SPACECRAFT WITH 


C ONTROL-MOMENT - GYROSC OPE 
STABILIZATION 

By Claude R. Keckler, Robert G. Kyle, Ralph W. Will, 
and Charles T. Woolley 
Langley Research Center 

SUMMARY 

A description of a real-time digital-computer-hardware simulation of an Apollo- 
Telescope-Mount (ATM) mission is presented. The simulation has been employed to 
evaluate the performance of a full-scale control-moment-gyroscope (CMG) system in 
providing stabilization for the ATM cluster. 

For this evaluation, a full-scale CMG system is linked to a real-time digital- 
computer program representing the ATM cluster. This program consists of the space- 
craft rigid- and flexible-body dynamics, internal and external disturbances, and the base- 
line CMG system momentum-feedback control law with isogonal correction. 

System-stability and limit-cycling problems were encountered during initial opera- 
tions. These problems are discussed and possible solutions are outlined. 

INTRODUCTION 

Theoretical research on the feasibility of using control-moment gyroscopes 
(CMG's) as the primary attitude -control system for long-term manned missions has indi- 
cated that CMG's offer significant advantages over other more conventional control 
concepts. 

The CMG is a double-gimbal device, which introduces torques by the precession of 
a constant-speed inertia wheel. A control system employing CMG's is capable of pro- 
viding the wide torque range necessary for performing all of the control tasks, such as 
maneuvering and fine pointing, required for long-term manned missions, such as the 
Apollo Telescope Mount (ATM). Furthermore, a CMG control system minimizes the 
optics -contaminant problems of reaction-jet systems, since CMG momentum dumping can 
be accomplished through gravity -gradient desaturation as shown by information obtained 
from H. F. Kennel at the Marshall Space Flight Center. Because of their unique 



characteristics, CMG systems also offer significant weight savings for long-term mis- 
sions over conventional control systems. 

To verify the results from preliminary CMG research, a real-time digital- 
computer simulation of the ATM vehicle was developed. In the course of the CMG 
developmental research, several analog- computer simulations have been developed and 
utilized effectively. The problem of simulating the complete ATM mission is much 
larger and more complex than any of the preceding simulations. Because of many alge- 
braic computations, a large dynamic range, and high resolution, digital computation was 
required. Hardware complexities and nonlinearities of the CMG system make a complete 
and realistic simulation of this hardware extremely difficult. Accurate representation of 
the CMG system characteristics is best provided by including the hardware in the simu- 
lation. This inclusion imposes the requirement of real-time operation. The first such 
real-time simulation employing a digital computer linked to a full-scale CMG control sys- 
tem has been developed and successfully utilized at the Langley Research Center for the 
Marshall Space Flight Center ATM mission. 

The general objectives of the static test program at the Langley Research Center 
were (1) to determine the CMG hardware characteristics, (2) to evaluate the hardware in 
the spacecraft control loop, and (3) to refine the response characteristics of the CMG 
system. 

This report presents a detailed description of the simulation used in the static test 
program in the form of flow charts, equations, hardware description, and discussions of 
simulation and system problems and solutions. 

SYMBOLS AND ABBREVIATIONS 

Symbols 

intermediate quantities defined by equations (34) 

3x3 direction-cosine input matrix to isogonal correction law 
modified direction cosine 

isogonal total CMG momentum component, lbf-ft-sec (N-m-sec) 
error of the isogonal distribution law 
error of the rotation law 
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crew-motion-disturbance force applied at ATM crew location, Ibf (newtons) 
crew-motion-distui’bance force applied at CSM crew location, lbf (newtons) 
crew-motion-disturbance force applied at S-IVB crew location, lbf (newtons) 
control-torque components of CMG system, lbf-ft (N-m) 

total generalized disturbance applied to flexible body, lbf-ft; lbf (N-m; newtons) 
individual CMG momentum magnitude, lbf-ft-sec (N-m-sec) 
command-momentum components, lbf-ft-sec (N-m-sec) 
momentum-error components, lbf-ft-sec (N-m-sec) 
momentum components of CMG system, lbf-ft-sec (N-m-sec) 

96 x 96 total-body generalized inertia-mass matrix, lbf-ft-sec2 (N-m-sec^) 
principal moments of inertia of body, slug-ft 2 (kg-m2) 
inertia, lbf-ft-sec 2 (N-m-sec2) 

96 x 96 total-body generalized stiffness matrix, lbf-ft/rad (N-m/rad) 
gain of gimbal-angle filter compensation, seconds 

actual and nominal isogonal distribution- law gains, rad/lbf-ft-sec 2 ® 

(rad/N-m-sec2) 

rotation-law gain for isogonal correction law, lbf-ft (N-m) 

effective inner-loop gain, rad/lbf-ft-sec 2 (rad/N-m-sec2) 

steering-law gains for outer and inner gimbal of CMG, rad/lbf-ft-sec 2 
(rad/N-m-sec^) 

proportional torque-command-attitude gains of CMG system, lbf-ft/rad 
(N-m/rad) 
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proportional torque-command-rate gains of CMG system, lbf -ft /rad/sec 
(N-m/rad/sec) 

distance of application of ATM crew-motion disturbance from cluster center 
of gravity, feet (meters) 

distance of application of CSM crew-motion disturbance from cluster center 
of gravity, feet (meters) 

distance of application of S-IVB crew-motion disturbance from cluster center 
of gravity, feet (meters) 

crew-motion-disturbance moment applied at ATM crew location, lbf-ft (N-m) 
filtered command moment component, lbf-ft (N-m) 

crew-motion-disturbance moment applied at CSM crew location, lbf-ft (N-m) 

crew-motion-disturbance moment applied at S-IVB crew location, lbf-ft 
(N-m) 

command-moment components of CMG system, lbf-ft (N-m) 
total external-disturbance moment components, lbf-ft (N-m) 
total internal-disturbance moment components, lbf-ft (N-m) 

96 generalized normal coordinates 
96 generalized normal rates, rad/sec 
96 generalized normal accelerations, rad/sec^ 

3x7 modal-participation matrix 
determinant of [e] 

intermediate quantities defined by equations (38) 
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total torque applied to rigid body, lbf-ft (N-m) 
time, seconds 

7 x 18 modal-distribution matrix 
spacecraft axes 
inertial reference axes 
96 spatial coordinates 
96 spatial accelerations, rad/sec^ 
orbit reference axis 
CMG outer-gimbal angles, degrees 
CMG outer-gimbal rates, deg/sec 

continuous CMG outer-gimbal-rate commands, deg/sec 
periodic isogonal correction for outer-gimbal-rate commands, deg/sec 
CMG inner-gimbal angles, degrees 
CMG inner-gimbal rates, deg/sec 

continuous CMG inner-gimbal- rate commands, deg/sec 

periodic isogonal correction for inner-gimbal-rate commands, deg/sec 

flexible-body attitude angles, radians 

inertial to orbital reference angle, radians (see fig. 6) 
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Subscripts: 

x,y,z 

1,2,3 


inertial to body reference angle, radians (see fig. 6) 
bending-filter damping coefficient 

angle between orbit plane and sunline, radians (see fig. 6) 
96 x 96 modal-participation matrix 
total-body modified Euler angles, radians 
commanded-body modified Euler angles, radians 
flexible-body modified Euler angles, radians 
rigid-body modified Euler angles, radians 
rigid-body modified Euler rates, rad/sec 
total -body angular rates, rad/sec 
commanded-body angular rates, rad/sec 
flexible-body angular rates, rad/sec 
orbit rate, rad/sec 
rigid-body angular rates, rad/sec 
rigid-body angular accelerations, rad/sec2 
96 x 96 modal-frequencies matrix, rad/sec 
bending-filter cutoff frequency, rad/sec 

components along spacecraft axes 
individual CMG in system 
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n integer 

T 

A [ j denotes a transpose matrix. 

Abbreviations 

ATM Apollo Telescope Mount 

ADC analog -to -digital converter 

CMG control-moment gyroscope 

CSM command and service module 

DAC digital-to-analog converter 

LM lunar module 

MDA multiple docking adapter 

MSFC George C. Marshall Space Flight Center 

S-IVB fourth stage of Saturn V launch vehicle 

SIMULATION CONCEPT 

Theoretical analyses of CMG systems have proven that these systems are capable 
of performing all the control tasks required by extended manned missions. (See ref. 1.) 
However, because of the inherent nonlinearities of CMG hardware as well as the complex- 
ity and inadequacy of CMG math models, these analyses have been limited to the assump- 
tions of idealized characteristics for the CMG control- system hardware. To evaluate the 
impact of the nonlinear characteristics of the hardware on the control capability of the 
systems, a simulation of a typical manned spacecraft application was developed and linked 
to a full-scale CMG system. The ATM cluster, shown in figure 1, was chosen as being 
representative of all near-future manned spacecraft. 

The ATM cluster consists of a repressurized S-IVB booster with a multiple docking 
adapter, an Apollo command and service module, and a modified lunar module. A gim- 
baled telescope experiment package is mounted in a rack attached to the LM. Three-axis 
spacecraft pointing to approximately 0.1° and 1 arc min/sec is required to support this 
experiment. 
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Analog simulations involving CMG systems have been evaluated but were found to 
be inadequate because of the large dynamic-range and resolution requirements associated 
with the ATM problem. Also, because of the complexity of the complete ATM simulation, 
a large number of algebraic computations were necessary which are not effectively per- 
formed on analog equipment. Therefore, it was apparent that a digital simulation was 
needed to meet these requirements; and since the simulation would be utilized in conjunc- 
tion with CMG hardware, it must operate in real time. Digital hardware, while capable 
of efficient algebraic computation and high accuracy, is not as efficient as analog hard- 
ware in performing integrations. It thus became necessary to find an efficient and fast 
integration algorithm to meet the accuracy requirements of the ATM problem and still be 
compatible with real-time operation. In addition, the iteration interval required to pro- 
vide a proper representation of the frequencies of the spacecraft and control system had 
to be established. A non-real-time digital program with a variable computing interval 
was developed and utilized to establish these parameters. 

Preliminary-response and accuracy analyses using the non-real-time program indi- 
cated that a fourth-order Runge-Kutta integration of the spacecraft and control-system 
equations with an iteration interval of 0.0312 second was required for proper frequency 
representations. Since the CMG hardware is restricted in frequency to 2.5 hertz 
(1 hertz = 1 cps), the 0.0312 second computing interval allows 13 samples in the period 
of the highest frequency content and is sufficient to accurately represent the system 
response. A typical iteration frame is shown in figure 2. As can be observed, only 
0.015 second is required to actually perform the computations; however, to meet real- 
time computer time-sharing requirements, a 0.0312-second interval is used. At the 
beginning of the computation interval, CMG gimbal rate and position information from the 
hardware are read into the program, and the integration step is initiated. At the end of 
each computation interval, the control-law commands are updated and transmitted to the 
CMG hardware. 

The ATM cluster is represented by seven primary modes of flexibility. Spacecraft 
rate and attitude output of the computer program represent motions of the ATM module 
which houses the telescope. All body-angular-rate and attitude sensors are assumed 
linear and their outputs are represented by the rate and attitude outputs of the digital 
program. The control law used in the simulation represents the baseline CMG flight sys- 
tem, that is, the momentum -feedback concept (ref. 2) with isogonal correction (developed 
by H. F. Kennel at MSFC). 

The simulation equations were programed in Fortran IV language on a Control Data 
series 6000 computer. This program is linked to the full-scale CMG system by trunk 
lines and by 14-bit plus sign-conversion equipment. (See fig. 3.) Simulation scaling is 

also shown in figure 3. 
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SIMULATION DESCRIPTION 


A block diagram of the simulation is presented in figure 4. The digital program 
includes control and disturbance torques, rigid-body dynamics, vehicle flexibility, and the 
CMG system control law. The simulation can be divided into two major parts: the outer 
or vehicle loop and the inner or momentum-feedback loop, which is indicated by heavy 
lines in figure 4. The control torque generated by the CMG's and external- and internal- 
disturbance torques are summed, and this result is used as an input to the spacecraft 
rigid-body equations. The output of this computation is the rigid-body contribution to the 
body-attitude and angular rate errors. In addition, the CMG control torques and internal 
disturbances are used to determine the rate and attitude errors caused by the spacecraft 
flexibility. The external-disturbance torque is considered to be distributed over the 
entire vehicle, and the frequency is low enough so as not to excite the flexibility modes 
of the spacecraft. Therefore, the flexibility contribution to the spacecraft errors 
resulting from the external torques is neglected. The sum of the flexibility- and rigid- 
body-rate and attitude errors is used as spacecraft-sensor-output signals, since these 
sensors were assumed to be linear. By using these output signals, a required or com- 
mand torque to correct for the spacecraft errors is generated. This command torque is 
passed through fourth-order bending filters to remove the higher frequency content. 

This procedure is necessary since the CMG system, because of its limited bandpass, 
would not respond properly to the high-frequency commands and would result in excessive 
wear of the CMG hardware. The filtered command torque is limited and integrated to 
yield the command momentum. This command momentum is summed with the CMG sys- 
tem momentum to yield the change in CMG system momentum required to correct for the 
spacecraft errors. The CMG gimbal rates required to provide the desired change in 
CMG momentum are then calculated by using the steering law. These commands are 
summed with isogonal corrections and transmitted to the hardware. The response of the 
hardware is monitored by the CMG gimbal sensors. The gimbal-sensor signals are then 
used to calculate the control torques generated by the hardware; thereby the simulation 
loop is closed. 

The governing equations for the simulation are presented herein in the order used 
by the digital computer during each iterative cycle. 

SPACECRAFT-TORQUE EQUATIONS 
Control-Torque Equations 

The control torques applied by the CMG system are computed from the following 

scalar equations. The CMG gimbal rate and position information required by these equa- 
tions are supplied by the CMG hardware gimbal sensors. The spacecraft-rate information 
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is obtained from the previous computer iteration or from the initial conditions since all 
vehicle sensors are idealized in this simulation. By using the MSFC coordinate system 
(see appendix), the CMG system control torques are expressed as 


Gy — H 
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These control-torque components are input to both the rigid-body equations and the 
flexible-body equations for the solution of total-body response. 

Disturbance -Torque Equations 

The simulation incorporates both the internal- and the external-disturbance torques 
to which the spacecraft is subjected. 

The internal disturbances consist of crew motions in the workshop, in the CSM, and 
in the LM and are represented in the form of three moments and three forces. The dis- 
turbances from crew motions employed in the simulation are either calculated from typi- 
cal crew-motion profiles, which were obtained by the Space Support Division of the Sperry 
Rand Corporation under contract to MSFC, programed on the computer or from actual 
ground-based crew-activity simulation data recorded on magnetic tape for input into the 
simulation. The equations representing the resultant torque from crew motion are 


m l v = M Ax + M Wx + Mj x + F Wy L Wz + F Jy Lj z - F Ay L Az 
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( 2 c) 


M t z = M A z + M W Z + M J Z + F A y L A x + F J y L J x “ F W y L W x 

where the subscripts x, y, and z refer to the body axes and the subscripts A, W, 
and J refer to the ATM, S-IVB workshop, and CSM areas, respectively. The L A , Lj, 
and terms are the moment arms from the crew-motion point of application to the 

cluster center of gravity and are always taken to be positive. (See fig. 5.) In addition, 
because of the cluster configuration, the distances along the Y-axis are relatively small; 
and thus, for the simulation, all crew-motion activity is limited to the X-Z plane. 

The forces and moments of equations (1) and (2) also act as disturbances (GE) n 
for the flexible body and are defined as 

(GE), = G x + M Ax ' 

(GE)2 = G y ♦ M Ay 
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(GE) 13 = Mj^ 
(GE) 24 = Mj y 

(GE) 15 = M Jz 
(GE) 16 = Fj x 
(GE)^ = Fj y 
(GE)ib = F Jz ^ 


applied at CSM location 


(3c) 


The external disturbances used in the simulation are limited to the gravity-gradient 
torques since these torques are the major orbital environment disturbance encountered 
by the vehicle for its particular mission profile. Therefore, from information obtained 
by H. F. Kennel at the MSFC, the torques imposed by external disturbances are 
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The inertial, orbital, and body-coordinate systems and their relationship used in devel- 
oping these equations are shown in figure 6. 


SPACECRAFT EQUATIONS OF MOTION 


Since the ATM cluster is a highly flexible structure, it becomes necessary to define 
both its rigid and flexible contributions to the total angular motions of the spacecraft. 

Rigid-Body Equations 

The equations of motion for the rigid-body contribution (assuming negligible prod- 
ucts of inertia) to the angular motion of the spacecraft can be expressed as follows: 
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(5c) 
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where the total-torque vector ^Tj^ is defined from equations (1), (2), and (4) as 

(T R ) = {G> + (Mj + (M e ) (6) 

These equations are solved for the rigid-body accelerations which are integrated to 
yield the rigid-body contribution to the vehicle angular rates. 

The rigid-body modified Euler angles 0 r, 0r, and 0 R are derived from the 
integration of the following differential equations: 
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Flexible-Body Equations 

The ATM vehicle is represented in this simulation by a modal model truncated to 
seven dominant modes. The contribution of the flexibility to the body-angular-rate and 
attitude errors is computed from the equations presented in this section. These equa- 
tions (see ref. 3) represent a model of the flexibility characteristics of the ATM cluster 
as obtained from the vibration analysis done under contract to MSFC during the system- 
design phase. 

The modal frequencies and normal modes for the vehicle are found by considering 
the following generalized equations of motion: 

[l]{x} + K {x} = {GE} (8) 

This set of equations can be uncoupled by assuming a solution of the form 


{*} - 03 <n> 


(9) 



where the [V] matrix represents the varying degree of participation of each normal 
coordinate in the flexible response. 

Since the participation matrix is constant, equation (9) may be differentiated twice 
to yield 



Inserting equations (9) and (10) into equation (8) forms 


HH(n} +[k]0{n} = {ge} 


This equation may be premultiplied by 



to become 


+H t [k][*]{n} =B t {ge} 


( 10 ) 


( 11 ) 


( 12 ) 


with the conditions 


and 






(13) 


(14) 


which the modal-participation matrix j~dTj must satisfy. Here, [}l{| is a diagonal unity 
matrix, and I s a diagonal matrix containing the squares of the modal frequencies. 

By considering structural damping, which for the ATM cluster is assumed to be 1 percent 
of critical damping, the set of equations for the normal modes becomes 


{n} + 2(0.01) (X) {n} 



{n} = (>] T {GE} 


(15) 


where the matrix distributes the generalized force terms over the set of modal 

equations. 


The flexible modes which were represented actively are those which fall within the 
range of the bending filter, as obtained from work done under contract to MSFC, and 
control-system bandwidths and are as follows: 
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Mode number 

Modal index 

Frequency, Hz 

1 

1 

0.149 

3 

2 

.161 

5 

3 

.414 

11 

4 

.540 

14 

5 

.638 

23 

6 

2.173 

25 

7 

2.52 


For the chosen representative flexible modes, equation (15) is partitioned to 

N-j + 0.02co i N i + wJn. = [u] . {GE} (16) 

J J J J J J 

where j = 1,2,. ..,7 is the modal index and jjj] is the 7 x 18 modal-distribution 
matrix, given in table I, which is a partition of the elements of the jjhj T matrix. 

Equation (10) can be expressed as 


% k ~ ^k,j N j < 17 ) 

j = l)2,...,7 

k=x,y,z 

or 

{%} =[p]{n} (18) 

where the Qp] matrix, presented in table II, is a 3 x 7 matrix comprised only of the 
columns of the (jQ matrix indicated in the summation. The flexible contribution to the 
body attitude is then 


{«} = .[{%}* 


(19) 


Total-Body Response 
The total -body rates are given by 


{d} = {%} + {fip} 


(20) 
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The flexible-body attitudes, which are assumed to be small with respect to those of 
the rigid body, are transformed into Eulerian coordinates by using the following approxi- 
mation in terms of the rigid-body Euler angles: 
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Finally the total-body modified Euler angles representative of spacecraft attitude 
are given by 
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CONTROL-SYSTEM EQUATIONS 
Momentum -Feedback Law 

The Euler angles and body angular rates developed in the equations of motion are 
used as the spacecraft signals in the outer-loop control law. The required torque M c 
for spacecraft control is expressed as 
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The commanded torque is filtered by fourth-order filters which prevent the high- 
frequency content caused by the spacecraft flexibility from being included in the CMG 
system gimbal-rate commands. This filtering of the command torque is required by the 
bandwidth limitations of the CMG system. The filter transfer functions used in the sim- 
ulation are the baseline flight-system filters and are given by 
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(24a) 



(24b) 


(24c) 


These resulting filtered moments are limited to ±200 lbf-ft (271.2 N-m) per axis 
and then integrated to yield the commanded momentum as shown in figure 7. The com- 
mand momenta H c constitute the inputs to the inner loop. It should be noted that H c 
has a limit of ±6000 Ibf-ft-sec (8136 N-m-sec). 

In the control-system inner loop, as shown in figure 4, the CMG momentum vector 
is calculated from the relationship 



(cos 0J COS O!^ 

- (cos 03 sin 

- (sin 0 2 ) 

sf 

II 

O 

(cos @2 cos 

- (cos 0 i sin 

“(sin 0 3 )j> 


(cos 03 cos a^\ 

- (cos 02 sin 0 * 2 ^ 

- (sin 0i) 


The required change in momentum desired from the CMG system to achieve space- 
craft control is then 


{H e } = {H c } - { H } 


(26) 
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The CMG gimbal rates to produce the required momentum change are computed by 
the steering law of reference 2. The steering-law transformation results in the following 


commanded CMG gimbal rates: 

0 C j = K SL^( - ^e x cos a \ sin 01 + H ey sin sin 0 J - H ez cos 0 ^ (27a) 

0c 2 = K SL^(“ H e x cos 02 " H e y cos 0 % sin 0 2 + H e z sin sin p 2 ) (27b) 

0c 3 = K SL^( H e x sin «3 sin 03 - H e y cos 03 " H ez cos a 3 sin 0 3 ) (27c) 

x = K SL a (" H e x sin ^ - H 6y cos (27d) 

“c 2 = K SL a (" H e y sin a 2 - H ez cos a 2 ) (27e) 

“c 3 = K SL a (‘ H e x cos 0:3 - Hg z sin Oj) (27f) 


where K SLq! and K SL ^ are constants with units of rad/lbf-ft-sec 2 (rad/N-m-sec2), 

and H e is limited to ±33.33 lbf-ft-sec (45.2 N-m-sec) per axis, as in the flight-system 
design. 


Isogonal-Correction Law 

Additional gimbal-rate commands to avoid antiparallel conditions are generated by 
the isogonal correction logic obtained from H. F. Kennel of MSFC. The isogonal com- 
putations are performed onboard the spacecraft by a digital computer and, as in the sim- 
ulation, are updated every 0.25 second. 

The isogonal correction consists of two parts: a distribution law which equally 
distributes the individual CMG momentum vectors about the instantaneous total-system 
momentum vector and a rotation law which minimizes the CMG inner-gimbal angles. 

The input to the digital computer consists of the direction cosines of the momentum 
components of each individual CMG about each spacecraft axis as follows: 


(28a) 

(28b) 
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E lx = cos cos aj 
Ei y = -cos 01 sin c^i 

E lz = -sin 0! 


(28c) 



E 2x = “Sin |3 2 

(28d) 

E 2 y = cos /3 2 cos a 2 

(28e) 

E 2z = -cos /3 2 sin 0^2 

(28f) 

e 3 x = -cos /3 3 sin 

(28g) 

E 3y = " Sin 03 

(28h) 

E 3z = cos 0 3 cos 

(28i) 


The distribution- and rotation-law gains used by the isogonal computer are K D p 
and K r , respectively. Additional direction cosines not available from the hardware are 
computed as follows: 



The total CMG system momentum Ep and the total components about each spacecraft 
axis are 


Erp — El + Eo + Eq 
A x A x °x 


E Ty = E ly + E 2 y + E 3y 
E T Z = E l z + e 2 z + e 3 z 


E™ = E™ + E™ + Em 
1 l x A y 1-2 


A 


Ep = ^Ep 


(30) 


and the solution of the determinant of direction cosines for the angular-momentum vec- 
tors becomes 
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q = E lx (E 2y E 3z - E 3y E 2z ) + E 2x (E 3y E lz - E ly E 3z ) + E 3 x (e % E 2z ■ E 2y E lz ) 

(sgn q = +1 for q § 0 
ysgn q = - 1 for q < 0 


( 31 ) 


where q = 0 defines a right-handed system and q < 0 defines a left-handed system of 
direction-cosine vectors. This information is used to determine the sign of the 
distribution-law gain. The magnitude of the distribution-law gain X is variable and is 
determined from the magnitude of the total-system momentum vector Ep with maximum 
gain occurring in the region of antiparallelism ^Ep ~ lj 


X = 0 


for 

E T £ 0.25 

X 5 = 2E'p 

- 0.5 

for 

0.25 < E T = 0.75 

A = +1 


for 

0.75 < E T £ 1.25 

X — +3 . 5 

-2E T 

for 

1.25 < E t ^ 1.65 

X s= 4-0.2 


for 

1.65 < E t 


(32) 


The actual distribution-law gain expressed as a function of X thus becomes 

K d = K Dp [(sgn q)x] 


(33) 


The following computations yield the CMG gimbal rates required to maintain the 
isogonal distribution and to minimize the CMG inner-gimbal angles. The distribution- 
law calculation involves taking the dot product of each individual CMG momentum vector 
and the total-system momentum vector and commanding gimbal rates which will equalize 
these dot products while maintaining constant system momentum. The rotation law com- 
mands CMG precession rates to minimize mean inner-gimbal angles. These required 
precession rates are then transformed into CMG gimbal coordinates to yield gimbal-rate 
commands as follows: 


(BE) X = E 2x E 3x + E 2y E 3y + E 2 z E 3z 
(BE) y = E 1 x E 3x + E ly E 3y + E 1 z E 3z \ 



( 34 ) 
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( 35 ) 



(BE) Z - (BE), 


Ki 


■D 


ej = [(BE) X - (BE) z Jk d \ 

% = [ (BE) y " < be >x]k d 


H, = e T Eo + e, Eg 
A x 1 z A y °x 

% = % E h * e h E h 

H lz. = e I z E 2z + e I y E 3z 

Ho = 6j Eo + 6 t Ei 
^x L x °x A z A x 

H9 = Gt Eo + Gt Ei 
*y A x ^y A z A y 

Ho = Gt Eo + Gt Ei 
*Z L X 6 Z A Z A Z 

Ho = Gt Ei + Gt Eo 

J x l y l x l x *x 

Ho = Gt Ei + Gt Eo 

6 y A y A y lx z y 

Ho = Gt E -j + Gt Eo 

6 z A v A z A x ‘•z 

J J 

E lz( E lz) 2 ( E ly E Tx - E lx E Ty) 

K r 

2 

E 2 X ( E 2 X ) ( E 2 z e T y - E 2y E T Z ) 
K r 

E 3y( E 3y) ( E 3 x E T z " E 3z E T x ) 


Kt 


( Gf ] 


1 + 6r 2 



J 


(36) 


(37) 
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R t = Hi + e r Enn 1 
A x l x r A x 

R ly = H ly + e r E Ty 

Ki •— Mi + G y> Jbj m 

x z i z A z 

r 2 x = h 2 x + e r E T x 

R 2y = H 2y + e r E Ty / 

Ro = Ho 4* G-pErp 

^z x z 

Ro ~~ Ho 4* GpEm 
°x °x X x 

Ro = Ho 4* 6-pErp 

6 y 6 y 1 A y 

Ro = Ho 4- 0Y» E rp 

°z °z x z 

J 

The CMG gimbal rates commanded by the isogonal computer are 

= - E i z ( E ly R l x • E lx R ly) 

^CI 2 = - E 2 x ( E 2 z R 2y - E 2y R 2 z ) 

“ci 2 = ( E 2 x ) 2E 2 x ( E 2y R 2y + E 2 z r 2 z ) " R 2 X 

^CI 3 = - E 3y( E 3 x R 3 z " E 3 z R 3 x ) 

"CI 3 = ( E 3y) 2E 3y( E 3 z R 3 z + E 3 x R 3 x ) - R 3y 


(38) 


(39) 


The isogonal-correction commands are limited to ±l°/ sec as in the flight system. 

The steering-law commands (eqs. (27)) are summed with the isogonal-correction 
commands (eqs. (39)) to yield total CMG gimbal-rate commands which are then fed to the 

hardware. The total CMG gimbal-rate commands are limited to a maximum of ±3.5°/sec 
to avoid overstressing the CMG rotor of the flight hardware. 
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CMG HARDWARE DESCRIPTION 


The hardware utilized in the real-time digital-hardware simulation of the ATM 

flight consists of three full-scale CMG's with a momentum capacity of 1000 lbf-ft-sec 
(1356 N-m-sec) per gyroscope. The CMG's are statically mounted on a unitized rigid 
test stand. (See fig. 8.) Each CMG is a double-gimbaled device which produces torque 
by precession of a constant-speed inertia wheel. Each gimbal is provided with a direct- 
current torquer acting through a compound planetary gear train and with gimbal-rate and 
position monitors. Excitation to the spin motors and gimbal torquers is provided through 
the electronics located in the master -control room. (See fig. 9.) Link between the CMG 
hardware and the digital computer is effected through the signal-buffer system in the 
master -control room. 

The CMG gimbals are driven and controlled by a rate servo loop (fig. 10) which 
uses a tachometer to monitor gimbal-torquer speed rather than actual gimbal rate as the 
feedback element. This servo-loop design was predicated on the need for a high-gain 
servo loop to provide adequate transient response and on the assumption that the gear 
train provided a stiff link between the torquer and the gimbal. 

From figure 11, it can be seen that the servo loop using the motor tachometer 
results in a bandwidth of approximately 150 rad/sec at the 90° phase point. This band- 
width would appear to be more than adequate to meet all the control requirements for 
manned vehicles such as the ATM. However, in monitoring the actual gimbal rate by 
means of a tachometer located after the gear train, large bandwidth discrepancies were 
noted. As can be observed from figure 12, the servo-loop bandwidth based on actual gim- 
bal rate is only about 20 rad/sec at the 90° phase point. The drastic reduction in fre- 
quency response was thought to be caused by excessive cross coupling between gimbals or 
by the fact that a rigid link between motor and gimbal did not exist as originally was 
assumed. Therefore, tests were made with one gimbal mechanically locked, thereby 
virtually eliminating the gimbal cross -coupling effects. The servo-loop bandwidth was 
not appreciably affected (fig. 13) ; and cross coupling was thus eliminated as the major 
cause for the loss in frequency response. Furthermore, the gimbal-tachometer output 
was monitored to observe gimbal response to sinusoidal and step inputs. For sinusoidal 
inputs, it was noted that a frequency was reached at which no torque or resultant motion 
was transmitted through the gear train to the gimbal. It should be recalled that the ampli- 
tude of the excursion or travel of a system such as a gear train will decrease as the fre- 
quency of the sinusoidal input is increased since less time elapses between command 
reversals. Thus, for commanded frequencies at which the excursion amplitude is equal 
to or less than the gear-train dead band, no gimbal motion will occur. The frequency 
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response of the CMG is therefore directly limited by this gear-train effect. Windup 
effects of the gear train caused by low stiffness introduce a time lag between input and 
output as well as overshoots in the output. Other servo loops have been tested on the 
hardware with no gain in bandwidth or time constant again because of gear-train limita- 
tions. Gear-train improvement is therefore necessary for better CMG hardware 
response. 


SIMULATION PROBLEMS 

The real-time simulation is conducted with the aid of a Control Data 6600 digital 
computer. However, the CMG hardware used in this simulation is an analog device. 
Thus, conversion equipment to transmit signals between the digital and analog realms is 
required. The 14-bit plus sign- conversion equipment utilized in this simulation has a 
resolution of approximately 1 part in 16 000 (214). However, because of the momentum- 
feedback control-law requirement for high resolution in the gimbal-position signals, this 
equipment was found to be inadequate, even though it is the best commercially available 
equipment. This lack in resolution was demonstrated by a flicker in the last bit of the 
analog-digital converter through which the CMG gimbal angles are passed. This flicker 
had an equivalent amplitude of 0.6 millivolt; and because of the simulation gains (fig. 3), 
the flicker resulted in a rate-command error of one-fifth of the full-scale command. An 
attempt at eliminating the problem by programing the simulation inner loop on an analog 
computer proved unsuccessful since the simulation gain remained unchanged and the 
analog-digital converters were still required. The solution to this problem was found in 
the use of a compensation-filter scheme employed to compensate for the detrimental 
gear-train effects and is discussed in the following section. 

SYSTEM PROBLEMS 

The gear trains furnished with the CMG hardware exhibited an output stiffness of 
47 000 lbf-ft/rad (63.7 kN-m/rad) and a backlash of over 6° at the input (0.1° at the out- 
put). This excessive dead band resulted in an unstable system at the system-design 
gains. As is evident from figure 14, a stable system was obtained for a maximum inner- 
loop gain (Ksl) of 0.000875 rad/lbf-ft-sec2 (0.000645 rad/N-m-sec2), which is about one 
order of magnitude lower than the ATM flight-design gain. (See ref. 2.) 

To meet the ATM mission pointing requirements, it became necessary to try to 
improve the CMG hardware. The corrective measures to improve CMG performance 
involved decreasing the gear -train backlash by plating the gears and increasing the gear- 
pivot-bearing preloads. A plot of average gimbal --rate limit cycling (with no spacecraft 

disturbances) for increasing inner-loop gain (fig. 15) illustrates the improvement. The 
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gear-train backlash is evidently still much too large for the system to perform all of the 
required control tasks. Thus, the resulting system performance was still unacceptable 
since the inner-loop gain was still below the design value and because limit-cycle ampli- 
tudes were still too large. (See fig. 16.) A compensation-filtering scheme by Willard W. 
Anderson of the Langley Research Center for the position signals was incorporated in the 
simulation. (See fig. 4.) This scheme maintains the actual gimbal position in the low- 
frequency region but alters the actual position on a high-frequency basis in calculating 
the system momentum. This alteration is accomplished by adding the commanded gimbal 
rate to the position signal and then filtering the sum. The position signal used by the 
computer is thus 

^ _ “gimbal + ^“command 
1 + K a s 


This compensation filter eliminates the high-frequency limit cycle induced by the phase 
lag associated with the actual gimbal position due to gear-train backlash. Thus, for the 
low-frequency region 


“gimbal + ^“command 
1 + K^s 


“gimbal 


and for the high-frequency region 


“gimbal + K a“command 
1 + KrS 


“command 


This scheme also allows for continuous updating of the gimbal position. Fig- 
ures 16, 17(a), and 17(b) present limit-cycle amplitudes for full inner -loop-design gains 
( K sl = 0.00525 rad/lbf-ft-sec2 (0.00387 rad/N -m-sec2)j for compensation gains K a 
of 0, 1.0, and 6.0 seconds. For values of K a > 6.0 seconds, the inner -loop -limit cycle 
is no longer noticeable. The benefits of the compensation -filter scheme are further 
evidenced by the data presented in figure 18. The compensation -filter gain was set at 
K a = 16.0 seconds and the system subjected to a three-axis maneuver of 0.05°. As can 
be seen, no noticeable inner -loop-limit cycle is present after completion of the maneuver. 
The apparent gimbal rate observed at the end of the maneuver in the data of figures 17 
and 18 is a result of recorder zero offset and not a product of the simulation. 

The difficulties encountered with the gear-train backlash in this simulation is a 
serious system problem primarily caused by the momentum-feedback-loop sensitivity to 
the small variations in CMG characteristics. 
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Subsequent simulations will include the control-moment gyroscopes mounted in 
torque-measurement fixtures to provide direct-control torque input to the computer. 

This combination will be mounted in a dynamic environment provided by a three-axis servo 
table (fig. 19) utilized to represent the spacecraft motions. Simulations to determine 
manual control feasibility will also be conducted by including a pilot control console 
(fig. 20) in the simulation loop. 


CONCLUDING REMARKS 

A description of a real-time digital-computer-hardware simulation of an Apollo- 
Telescope-Mount mission has been presented. Problems encountered during the investi- 
gation consisted of both simulation and system problems. 

The simulation problem was caused by the resolution limitation of the conversion 
equipment. This limitation introduced an erroneous control-moment-gyroscope (CMG) 
gimbal-rate command equal to one-fifth of the full-scale gimbal-rate command. 

The system problem was an excessive limit cycle of the spacecraft during the fine- 
pointing mode. This limit cycle was directly attributed to the CMG gimbal actuators. 
These actuators exhibited high-backlash and low-stiffness characteristics which resulted 
in an unstable system for the system design gains. 

The solution to both of these problems was found in the use of a compensation- 
filtering scheme for the position-feedback signals from the hardware. This scheme 
maintains the actual gimbal position in the low-frequency region but alters the actual 
position on a high-frequency basis in the calculation of the system momentum. The 
compensation-filtering scheme utilized in this simulation was effective in alleviating the 
sensitivity of the momentum-feedback control law to small variations in CMG character- 
istics. The employment of this scheme resulted in a stable system at the system design 
gains and thus permitted a complete system evaluation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 11, 1970. 
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APPENDIX 


COORDINATE SYSTEM 

The axis placement and positive gimbal-angle notation for the control-moment 
gyroscope oriented as suggested by researchers at the Marshall Space Flight Center are 
defined by the following sketch: 

Z 



Outer -gimbal angles are denoted with the symbol a and inner-gimbal angles with jS. 
For nominal inner-gimbal angles, /3j = 02 = j3 3 = 0°, the outer-gimbal rate is about 
the -Z spacecraft axis, about the -X spacecraft axis, and d 3 about the -Y 
spacecraft axis. 
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TABLE L- MODAL-DISTRIBUTION MATRIX 


4.40 x 10“ 4 

4.80 x 10" 6 

2.28 x 10- 5 

3.00 x 10 -5 
4.22 x 10" 3 

-9.00 x 10 -5 

7.28 x 10“ 4 

4.80 x 10" 6 
2.16 x 10- 5 

-2.00 x lO" 5 
-1.55 x lO" 3 
-2.60 x lO" 4 

7.28 x lO" 4 

4.80 x 10“ 6 

2.28 x 10" 5 
-2.00 x lO" 5 
-3.08 x 10" 3 

6.00 x 10“ 5 


-4.27 x 10“ 4 
-8.40 x 10" 6 
-2.28 x lO" 5 
-6.00 x 10" 5 
-4.41 x 10~ 3 
-8.90 x lO" 4 
-7.92 x 10 -4 
-7.20 x 10“ 6 
-2.16 x lO" 3 

4.00 x lO" 5 
1.57 x 10" 3 
-6.10 x lO" 4 
-7.90 x lO- 4 
-8.40 x 10" 6 
-2.40 x lO" 5 

4.00 X lO" 5 

3.11 x lO" 3 
-1.14 x lO" 3 


6.91 x 10" 3 
-7.44 x lO" 5 

2.41 x lO" 4 
-4.60 X lO" 4 
1.73 x lO" 3 

7.90 x lO" 4 
-6.20 x 10“ 3 
-3.72 x 10-5 
1.96 x lO" 4 
2.40 x lO" 4 
-9.62 x 10~ 3 
2.16 x 10- 3 
-6.11 x 10" 3 
-3.84 x lO" 5 
2.54 x lO" 4 
2.40 X lO" 4 
-2.47 x lO" 2 
-4.10 x lO" 4 


-2.21 x lO" 3 

1.50 x lO" 4 
-6.72 x lO" 5 

6.90 x lO" 4 

4.64 x 10~ 3 
-5.70 x lO- 4 

3.38 X 10” 3 

2.64 x lO" 5 
-5.04 x lO" 5 
-4.80 x lO- 4 

1.90 x lO" 3 
-1.57 x lO" 3 

3.37 x 10“ 3 

3.12 x 10-5 
-7.08 X lO" 5 
-4.80 x lO" 4 

5.98 x 10" 3 

4.00 x lO- 4 


-1.21 X lO" 4 
-8.44 x lO" 3 
-3.60 x lO" 6 
-2.99 x lO" 2 
5.20 x lO" 4 
6.07 x lO" 3 

2.44 x lO" 4 

6.01 x lO" 4 

-4.80 x lO" 6 

2.50 x lO" 2 
-2.00 x lO” 5 
-1.32 x lO- 2 

2.45 x lO" 4 
3.98 x lO" 4 

-7.20 X lO" 6 

2.50 x lO" 2 

3.90 x 10-4 
1.83 x 10-2 


2.21 x lO" 3 
-8.88 x lO" 5 
-1.89 x 10- 3 
7.00 x lO" 5 
-2.15 x lO" 2 
-4.28 x lO" 3 
2.57 x lO" 4 
-1.70 x lO- 4 

1.65 x lO" 3 
-6.20 x lO" 4 
-1.21 x lO" 2 
-1.29 x lO" 3 

1.13 x lO" 4 

4.80 x lO" 4 
-3.25 X lO- 3 
-6.50 x lO -4 
5.47 x 10-2 
9.55 x lO” 3 


-2.92 x lO -3 
-6.60 x 10~5 
-1.01 x lO" 2 

2.80 x lO" 4 

3.65 x lO -2 
-1.54 x 10-3 
-3.52 x lO" 4 
-6.48 x lO" 5 
-1.88 x 10“3 
2.30 x lO" 4 
2.25 x lO" 2 
-3.00 x lO -5 
-6.96 x IQ-5 

1.37 x 10-4 

5.41 x lO -3 
2.70 x lO -4 

-9.20 x lO -2 

2.38 x ID -3 



CO 

o 


TABLE II.- MODAL- PARTICIPATION MATRIX 



3.67 x 10-5 

-3.56 X 10-5 

5.76 X 10-4 

-1.841 x 10-4 

-1.010 X 10-5 

1.838 X 10-4 

-2.432 x 10-4 

r~i 

lHj 

ii 

4.00 x 10~ 7 

-7.00 X 10- 7 

-6.20 x 10-6 

1.250 x 10-5 

-7.031 x 10-4 

-7.400 x 10-6 

-5.500 X 10-6 


1.90 x 10“ 6 

-1.90 X 10“ 6 

2.01 x IQ" 5 

-5.600 X 10-6 

-3.000 x 10- 7 

-1.573 x 10-4 

-8.388 x 10-4 







31 








Attenuation ~ 1/80 



Figure 3.- Simulation concept and scaling. 
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Figure 8.- CMG static test setup. 
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Figure 9.- Master-control room for CMG system. 
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Figure 11.- CMG frequency response from motor tachometer. 
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Figure 12.- CMG frequency response from gimba! tachometer. 
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Figure 13.- CMG frequency response from gimbal tachometer with no cross coupling. 
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